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Abstract 

The genus Sceloporus is one of the largest genus of lizards in North and Central America, with 22 species groups. 
Among these, the torquatus group has a notably wide geographic distribution with populations occurring from southern 
United States to Guatemala. In spite of the taxonomical work done with the group, some problems remain unsolved. We 
therefore obtained the phylogeny of the torquatus group, based on 925 bp of the ribosomal 16S gene, 912 bp of the ribo- 
somal 12S gene, and 893 bp of the ND4 gene, for a total of 54 specimens of 25 taxa. The genes were analyzed, both sep¬ 
arately and combined, by means of maximum parsimony and Bayesian inference analyses. The subspecies of .S', serrifer 
did not form a monophyletic group. The sequence data refuted the morphological evidence that suggested that .S', s. 
plioporus and S. cyanogenys are closely related to S. s. serrifer and to S. s. prezygus. Regardless, these last two were 
recovered as sister taxa. Moreover, evidence was found that S. ornatus does not form a monophyletic group, and that S. 
ornatus ornatus and S. oberon are a single species, despite their marked differences in coloration and scutelation. In addi¬ 
tion, the non-monophyly of S. mucronatus was confirmed and the phylogenetic relationships of its different species were 
determined. At the same time, the subspecies of S. dugesii were recovered as a monophyletic group, refuting the non- 
monophyly of this taxon suggested in the phylogenetic hypothesis of the entire genus. 
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Resumen 

El genero Sceloporus es uno de los generos mas grande de lacertilios de Norte y Centroamerica, con 22 grupos de espe- 
cies. Entre estos, el grupo torquatus tiene una amplia distribution geografica, con poblaciones que ocurren desde el sur 
de Estados Unidos hasta Guatemala. No obstante los trabajos taxonomicos realizados hasta ahora con el grupo, algunos 
problemas permanecen sin resolver. Por esa razon, obtuvimos la filogenia del grupo torquatus, basados en 925 pb del gen 
ribosomal 16S, 912 pb del gen ribosomal 12S y 893 pb del gen ND4, para un total de 54 especimenes de 25 taxa. Los 
grupos de datos fueron analizados separadamente y en conjunto, por medio de maxima parsimonia e inferencia bayesi- 
ana. Las subspecies de S. serrifer no fueron recuperadas formando un grupo monofiletico, los datos refutan la evidencia 
morfologica que sugiere que S. s. plioporus y S. cyanogenys se realcionan con S. s. serrifer y con S. s. prezygus, sin 
embargo estos dos ultimos si son recuperados como taxa hermanos. Asimismo, encontramos evidencia que sugiere que 
las subspecies de S. ornatus no forman un grupo monofiletico, y que S. ornatus ornatus y S. oberon forman parte de una 
sola especie, a pesar de sus marcadas diferencias en coloration y escutelacion. Tambien, se confirmo la no monofilia de 
S. mucronatus y se determinaron las relaciones filogeneticas de sus distintas especies. A1 mismo tiempo, las subspecies 
de S. dugesii se recuperaron como un grupo monofiletico, lo cual refuta la no monofilia de este taxon como se habia sido 
sugerido en la filogenia previa del genero. 
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Introduction 


The genus Sceloporus (Squamata: Phrynosomatidae) is probably the best represented genus of small lizards in 
North and Central America, with ca. 70 species distributed from the Northern United States to Panama. 
Within the genus, 22 groups are distinguished, with approximately 80 species, of which some 30 species are 
viviparous (Sites et al. 1992; Wiens and Reeder 1997). 

The torquatus group (sensu Smith 1938) contains viviparous species widely distributed from the southern 
United States southward into Guatemala. Most of these species occur in mountainous areas with temperate 
environments. A few species are distributed in lowlands with tropical or semi-desert environments (Smith 
1936; Smith 1939; Sites et al. 1992). This group was proposed by Smith (1938), originally with the name of 
torquatus group, but later Smith (1939) proposed that the group should be called the poinsetti group, because 
the name of torquatus was at one time a secondary homonymy and under the nomenclatural rules of that time, 
the suppression of any homonym was considered permanent (Bell et al. 2003). Nevertheless different rules 
were proposed and Smith and Taylor (1950) revived torquatus as a group name, but this change was made 
official until 1961, when the new Code appeared. 

The torquatus group can be diagnosed by a series of characters, described mainly by Smith (1938; 1939), 
with exception of some characters described later by Wiens and Reeder (1997): wide separation between 
xiphisternal ribs (Wiens and Reeder 1997); no contact between frontal and interparietal scales; no contact 
between frontal and median frontonasal scales; median parietal scale present; a lip below tip of scales (Wiens 
and Reeder 1997); granular skin between scales (Wiens and Reeder 1997); dorsal, ventral, and lateral scales 
distinctly differing in size; lateral nuchal scales not well differentiated from dorsal nuchal scales; dorsal scales 
subequal in size; a distinct dark, light bordered nuchal collar; male belly patch with incomplete dark margin 
(Wiens and Reeder 1997); female belly patches absents; karyological characters and DNA sequences (Wiens 
and Reeder 1997). 

Since the proposal of the torquatus group by Smith (1938), the description of many new species and sub¬ 
species continued with traditional morphological characters. Nevertheless, when Wiens and Reeder (1997) 
proposed the Sceloporus phylogeny based on molecular and morphological characters, they detected a 
strongly supported conflict between DNA and morphological data. The study of Wiens et al. (1999), in which 
molecular data for the majority of the subspecies of S.jarrovii Cope in Yarrow were obtained, found that these 
subspecies were not monophyletic, and suggested a series of nomenclatural changes within which five evolu¬ 
tionary species were proposed. Wiens et al. (1999) also proposed that the divergence in coloration is possibly 
the result of sexual selection and habitat features. In other study carried out by Wiens and Penkrot (2002), 
concerning the delimitation of species using DNA and morphological characters, they used species of the 
torquatus group for the exemplification of a new protocol, and again a discordance between morphology and 
DNA was found, where two mtDNA clades were recognized as species that lacked diagnostic morphological 
characters. They established that in the torquatus group there is a particular pattern of morphological varia¬ 
tion, in which between-species differentiation is small relative to within-species, the worst combination for 
morphology-based delimitation. For this reason, the taxonomic status of some species of the torquatus group 
remains uncertain, taking into account that in the studies of Wiens and Reeder (1997) and Wiens et al. (1999), 
no DNA sequence data were available for many taxa. 

We undertook a phylogenetic investigation of the torquatus group using mtDNA sequence data from all 
taxa. The results of this study are reported herein. 


Materials and methods 

Taxon sampling 

Samples of liver and muscle tissue were obtained from the 17 taxa not included in previous studies on the 


54 ■ Zootaxa 1609 © 2007 Magnolia Press 


MARTINEZ-MENDEZ & CRUZ 



torquatus group, and for which no sequence data existed in DNA sequence data banks. Additional samples 
included a new record (5. sp. 1), a population of S. bulleri Boulanger (5. bulleri 2 + 5. bulleri 3), additional 5. 
cyanogenys Cope (5. cyanogenys 2), 5. minor Cope and a recently discovered population of 5. torquatus mel- 
anogaster Cope (Hemandez-Gallegos et al. 2003) (Fig. 1, Tab. 1). We included sequences generated by Wiens 
and Reeder (1997) and Wiens et al. (1999), which were obtained from GenBank (accession numbers in Table 
!)• 



FIGURE 1. Distribution of Sceloporus torquatus species -group in Mexico, south of United States of America and Gua¬ 
temala, based on Smith 1938, Wiens et al. (1999) and museum data. Dots represents localities sampled for this study and 
those reported for each specimens included from GenBank. Numbers represents the taxa included in the analyses: 1. S. 
bulleri; 2. S. cyanogenys', 3. S. cyanostictus', 4. S. dugesii dugesii; 5. S. d. intermedius; 6. S. insignis ; 7. S.jarrovii; 8. 5. 
lineolateralis', 9. S. macdougallr, 10. S. minor, 11. S. mucronatus aureolus', 12. S. mucronatus mucronatus; 13. S. mucr- 
onatus omiltemanus', 14. S. oberon; 15. S. ornatus caeruleus; 16. S. ornatus ornatus; 17. S. poinsettii; 18. S. serrifer 
plioporus; 19. S. serriferprezygus; 20. S. serrifer serrifer, 21. S. sugillatus; 22. S. torquatus binocularis; 23. S. torquatus 
melanogaster; 24. S. torquatus torquatus', Sceloporus sp. 1; 26. Sceloporus sp. 2. The abbreviations means: In United 
States of America: AZ.= Arizona, NM.= New Mexico, TX.= Texas; In Mexico: CHIS.= Chiapas, COAH.= Coahuila, 
NL.= Nuevo Leon, TMPS.= Tamaulipas VER.= Veracruz and YUC.= Ycatan. 

DNA isolation, PCR amplification and sequencing 

MtDNA was isolated from small quantities of liver and muscle (approx. 100 mg) following Fetzer’s 
(1996) extraction protocol with ammonium acetate. The target genes were amplified using the polymerase 
chain reaction (PCR; Saiki et al. 1998). The amplified regions correspond to a fragment of approximately 912 
base pairs of the ribosomal 12S (rRNA) gene using the primers tPhe and 12e (Wiens et al. 1999); 925 base 
pairs of the 16S (rRNA) gene using the primers 16SaR-L and 16Sd-Fl (Reeder 1995) and 893 base pairs of the 
ND4 gene that additionally included the complete portions of t-RNA-Flis, t-RNA-Ser and a portion of the t- 
RNA-Leu, amplified with the primers ND4 and LEU (Arevalo et al. 1994; Forstner et al. 1995). These genes 
yielded good results in other studies of Sceloporus (Benabib et al. 1997; Wiens and Reeder 1997; and Wiens 
et al. 1999). 
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TABLE 1. Species, localities, voucher specimen number and GanBank accession numbers for specimens evaluated in 
the Sceloporus torquatus species-group. The acronyms follow the nomenclature of Leviton et al. (1985) except for 
MZFC, which corresponds to the Museo de Zoologia of the Facultad de Ciencias of the Universidad Nacional Autonoma 
de Mexico (UNAM); MX, MZFC frozen collection, and IBFl, which corresponds to the Coleccion Nacional de Anfibios 
y Reptiles of the Instituto de Biologia of the UNAM. The numbers and letters after .S', oberon and S. minor correspond to 
the population and organism code with which they are identified in the study of Wiens et at. (1999). 


Species 

Locality 

Voucher 

GenBank accesion no. 





12S 

16S 

ND4 

Sceloporus bulleri 1 

Mexico: Jalisco: 1.0 km S Mas- 

cota 

IBH 18034 

DQ525887 

DQ525904 

DQ525865 

Sceloporus bulleri 2 

Mexico: Jalisco 

MX15-63 


EF608027 

EF608022 

Sceloporus bulleri 3 

Mexico: Jalisco 

MX 15-64 


EF608028 

EF608023 

Sceloporus cyanogenys 

United States: Texas: 

McMullen 

LSUMZ 48852 

AF15414 

AF000876 

AF154193 

Sceloporus cyanogenys 2 

Mexico: Nuevo Leon: 

Escobedo: 25.3 km NW 
Monterrey 

IBH 18051 

DQ525893 

DQ525910 

DQ525868 

Sceloporus cyanostictus la 

Mexico: Coahuila: 23.6 km S 

Monclova 

CM 147644 

AF154146 


AF154194 

Sceloporus cyanostictus lb 

Mexico: Coahuila: 1.0 km S 

San Lorenzo 

MZFC 7411b 

AF000825 

AF000865 

AF154195 

Sceloporus dugesii dugesii 

Mexico: Jalisco: Tapalpa 

UTA-R 23955 

AF154170 

AF000877 

AF154190 

Sceloporus duguessi interme- 
dius 

Mexico: Guanajuato: 2.0 km E 
Moroleon 

IBH 18002 

DQ525886 

DQ525903 

DQ525878 

Sceloporus duguessi interme¬ 
dins 2 

Mexico: Guanajuato: 2.0 km E 
Moroleon 

IBH 18004 

DQ525889 

DQ525906 

DQ525866 

Sceloporus insignis 

Mexico: Michoacan 

no voucher 

AF000806 

AF000846 


Sceloporus jarrovii 11a 

Mexico: Zacatecas: 24 km W 

Fresnillo 

CM 147650 

AF15173 


AF154209 

Sceloporus jarrovii lib 

Mexico: Zacatecas: 24 km W 

Fresnilo 

CM 147651 

AF15418 


AF154210 

Sceloporus jarroviii 10 

United States: Arizona: 

Cochise Co., near Portal 

LSUMZ 48786 

AF154163 

AF000881 

AF154208 

Sceloporus lineolateralis 

Mexico: Durango: nearPedri- 

cena 

MZFC 6650 

AF000807 

AF000847 

AF154211 

Sceloporus macdougalli 

Mexico: Oaxaca: Rincon 
Bamba, 35. 2 km SW Tehua¬ 
ntepec 

MZFC 7017 

AF000809 

AF000849 


Sceloporus minor 

Mexico: Tamaulipas: 17.7 km 
SW Ciudad Victoria 

IBH 18012 

DQ525891 

DQ525908 

DQ525872 

Sceloporus minor 13a 

Mexico: Zacatecas: 4.0 km W 
Concepcion del Oro. 

MZFC 10703 

AF154185 


AF154222 

Sceloporus minor 14a 

Mexico: San Luis Potosi: Colo- 
nia Insurgentes, 2.5 km W San 
Luis Potosi 

CM 147653 

AF154174 


AF154218 

Sceloporus minor 15a 

Mexico: San Luis Potosi: 14.1 

km E Ciudad del Maiz 

CM 147630 

AF154136 


AF154213 





to be continued. 
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TABLE 1. (continued) 


Species 

Locality 

Voucher 

GenBank accesion no. 

12S 16S 

ND4 

Scelopoms minor 17 

Mexico: San Luis Potosi: 22.8 

km E Matehuala 

CM 147679 

AF154148 


AF154231 

Sceloponis minor 3 

Mexico: Queretaro: 4.9 km S 
Ezequiel Montes 

MZFC 10736 

AF154138 


AF154198 

Scelopoms minor 4b 

Mexico: Queretaro: 1.0 km S 
Cadereyta 

MZFC 10738 

AF154175 


AF154200 

Sceloponis minor 5 

Mexico: Hidalgo: Barranca de 
los Marmoles W of Jacala 

CM 147625 

AF154142 


AF154201 

Scelopoms minor 6a 

Mexico: Hidalgo: Puerto de la 
Zorra, between Cuesta Colo- 
rada and Jacala on Hwy 85 

CM 147628 

AF154143 


AF154202 

Scelopoms minor 8 

Mexico: Tamaulipas: 16.9 km 
W Ciudad Victoria 

MZFC 10666 

AF154155 


AF154204 

Scelopoms mucronatus aure- 
olus 

Mexico: Oaxaca: Temazulapan 

IBH 18022 

DQ525884 

DQ525901 

DQ525875 

Scelopoms mucronatus 
mucronatus 

Mexico: Estado de Mexico: 
Ajusco Volcano: Ejido Capulin 

IBH 18008 

DQ525885 

DQ525902 

DQ525864 

Scelopoms mucronatus 
omiltemanus 

Mexico: Guerrero: Omiltemi 

National Park 

UTA-R 24004 

L41419 

L41469 

AF154233 

Scelopoms oberon 21b 

Mexico: Nuevo Leon: 9.0 km E 

San Roberto 

MZFC 8032 

AF000826 

AF000866 

AF154212 

Scelopoms oberon 24a 

Mexico: Nuevo Leon: 2.1 km S 
Santa Clara de Cienega 

CM 147675 

AF154157 


AF154228 

Scelopoms oberon 27b 

Mexico: Coahuila: N of El Dia¬ 
mante 

CM 147674 

AF154183 


AF154239 

Scelopoms oberon 28a 

Mexico: Coahuila: 22.3 km E 

San Antonio de las Alazanas 

CM 147641 

AF154160 


AF154234 

Scelopoms oberon 29a 

Mexico: Nuevo Leon: 2.5 km E 
San Isidro, turnoff for Laguna 
Sanchez 

MZFC 10698 

AF154147 


AF154236 

Scelopoms ornatus caeruleus 

Mexico: Coahuila 

JAM 652 

AF000814 

AF000854 

AF154240 

Scelopoms ornatus ornatus 

Mexico: Coahuila: Ojo Cali- 
ente N of Ramos Arizpe 

IBH 18041 

DQ525879 

DQ525896 

DQ525862 

Scelopoms poinsettii 

United States: Texas: Val Verde 

Co. 

LSUMZ 48847 

AF154176 

AF000883 

AF154241 

Scelopoms serrifer plioporus 

Mexico: Tamaulipas: Padilla, 
4.5 km NW Ciudad Victoria 

IBH 18014 

DQ525882 

DQ525899 

DQ525873 

Scelopoms serrifer plioporus 

2 

Mexico: Tamaulipas: Padilla, 
4.5 km NW Ciudad Victoria 

IBH 18015 

DQ525883 

DQ525900 

DQ525874 

Scelopoms serrifer prezygus 

Mexico: Chiapas: 2.5 km NW 
Teopisca 

IBH 18027 

DQ525880 

DQ525897 

DQ525870 

Scelopoms serrifer prezygus 

2 

Mexico: Chiapas: Ixtapa, 26.3 
km E San Cristobal de las 

Casas on Hwy 190 

IBH 18095 

DQ525881 

DQ525898 

DQ525875 





to be continued. 
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TABLE 1. (continued) 


Species 

Locality 

Voucher 

GenBank accesion no. 

12S 16S 

ND4 

Sceloporus serrifer serrifer 

Mexico: Yucatan: 13 km N 

Merida 

IBH 18020 

DQ525894 

DQ525911 

DQ525876 

Sceloporus serrifer serrifer 2 

Mexico: Yucatan: 20 km N Tiz- 

imin 

IBH 18123 

DQ525895 

DQ525912 

DQ525877 

Sceloporus sugillatus 30a 

Mexico: Morelos: Lagunas de 
Zempoala, W ofHuitzilac 

CM 147623a 

AF154187 


AF154242 

Sceloporus torquatus binocu- 
laris 

Mexico: Nuevo Leon 

MZFC 8033 

AF000827 

AF000867 


Sceloporus torquatus melan- 
ogaster 1 

Mexico: N of Estado de Mex¬ 
ico 

UTA-R 24016 

AF154179 

AF000890 

AF154244 

Sceloporus torquatus melan- 
ogaster 2 

Mexico: Estado de Mexico: 2.0 

km N Polotitlan 

IBH 18006 

DQ525892 

DQ525909 

DQ525863 

Sceloporus sp. 1 

Mexico: Nayarit: Sierra de 
Alica: Carretera Huajimin- 
Tepic 

MX 14-4 

EF608018 

EF608026 

EF608021 

Sceloporus sp. 1 

Mexico: Jalisco: Bolanos 

MX 13-24 

EF608016 

EF608024 

EF608019 

Sceloporus sp. 1 

Mexico: Jalisco: Carretera 
Bolanos-Tuxpan de Bolanos 

MX 13-80 

EF608017 

EF608025 

EF608020 

Sceloporus heterolepis 

Mexico: Jalisco: Cumbre de los 

Arrastrados 

IBH 18138 

DQ525890 

DQ525907 

DQ525869 

Sceloporus grammicus 

Mexico: Oaxaca: Sierra de 

Juarez 

UTA-R 23970 

L40457 

L41464 

AF154188 


PCR reactions in a Perkin-Elmer 2400 thermocycler had a final volume of 50 pi. The conditions for the 
PCR reaction for the different genes were: 12S with 45 cycles of 94° C for 30 sec, 53° C for 30 sec, and 12° C 
for 2 min; 16s with 45 cycles of 94° C for 30 sec, 50° C for 45 sec, and 12° C for 30 sec; and for the ND4, 35 
cycles of 94° C for 1 min, 50° C for 1 min, and 12° C for 1 min were performed. The first cycle of each of the 
amplification reactions included a denaturalization cycle of 94° C for 3 min, and the last cycle was completed 
with a cycle for final extension (two in the case of ND4) of 12° C for 5 min. The PCR products were purified 
using the QIA quick purification kit, and the resulting samples were sequenced by means of the automated 
sequencing service of the UNAM’s Instituto de Biologia, utilizing an ABI PRISM, 310 Genetic Analyzer 
(Applied Biosystems) automated sequencer. All DNA sequences obtained were deposited in the GenBank 
(Accession Nos. DQ525862-DQ525912) and are listed in Table 1. 

Sequence alignment 

The sequences obtained were compiled and edited in ProSeq 2.91 (Filatov 2002). Sequence alignment 
was carried out separately for each region, employing Clustal X (Thompson et al. 1997), using the default 
parameters, and, later, manually refined using the secondary structural models for the 12S and 16S (Orti and 
Meyer 1987). Consequently, alignment was again made using Clustal X, with apertures and gap extensions of 
15:6, 10:5, 6:3 and 3:1. The sequence regions, whose homologies by the nucleotide position were at a vari¬ 
ance, in differing penalizations, were considered ambiguous and were not included in the phylogenetic analy¬ 
ses. The ND4 protein-coding genes lacked insertions and deletions (indels) and were aligned by eye. Later, 
this codifying region was transferred to amino acids in order to check whether stop codons existed that could 
indicate the presence of pseudogenes. The tRNAs region adjacent to the ND4 was also aligned by eye, and 
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was later reanalyzed with Clustal X using different gap costs. For the total evidence analysis, the matrices 
aligned for each region were combined into a new matrix. Different number of terminals for each region were 
coded as missing data. 

Phylogenetic analysis 

The phylogenetic analysis of the molecular data used total evidence (Kluge and Wolf 1993). In order to 
detect possible areas of significant incongruence, the genes were also analyzed independently (Wiens 1998a). 

We did not test for the presence of phylogenetic signal. The signal is additive across different matrices and 
can dominate in a combined analysis in cases where the separate matrices have a very weak signal (Barrett et 
al. 1991; Wenzel and Siddall 1999). Because some of the sequences were obtained from GenBank, the sepa¬ 
rate matrices did not include the same taxa. Nevertheless, all the sequences were included in the combined 
analysis to maximize sampling. 

Maximum parsimony analyses (MP) were conducted in PAUP version 4.0b 10 (Swofford 1998) for the 
separate and the total evidence data sets. We used a heuristic search with tree bisection and reconnection 
(TBR) branch swapping and 1000 random sequence addition replicates. Characters were treated as unordered 
and equally weighted, and gaps were coded as missing data. Branches were collapsed if the maximum length 
was zero. Clade support was evaluated using nonparametric bootstrap proportions (BSP, Felsenstein 1985) 
with 1000 pseudoreplicates. BSPs proportions of <70% were considered to indicate poor support (Brandley 
and De Queiroz 2004). BSPs of =95% were interpreted as representing very strong support and from 70% to 
94% moderate support. 

Modeltest (version 3.07, Posada and Crandall 1998) was used to infer the best-fit model of evolution for 
the Bayesian inference (BI) analyses for each partition based on the Bayesian Information Criterion (BIC) 
method. The Hierarchical Ratio Test (hLRTs), although being the most popular method, is not the optimum 
strategy for choosing substitution models for phylogenies (Sanderson and Kim 2000; Posada and Buckley 
2004). 

BI analyses (Larget and Simon 1999; Lewis 2001) were performed for 12S, 16S, and ND4 matrices (with 
four partitions: codons + tRNAs) and for the combined data set with the six previous partitions, using 
MrBayes 3.0 b4 (Huelsenbeck and Ronquist 2001). Because MrBayes is limited to models with one, two or 
six base-substitution rate matrices, we used the GTR+I+G model for 12S, 16S and ND4 second and third 
codon positions instead of TrN+I+G model (best model obtained by Modeltest) because TrN has three param¬ 
eters. For tRNAs we used the GTR+G model instead of K81uf+G model, because K81uf also has three param¬ 
eters. For the first position codon in ND4 we used the HYK+G model inferred by Modeltest. In each analysis, 
four Markov chains were run, beginning with a random tree. The analysis used 2.0 x 10 6 generations with sam¬ 
pling every 1000 generations. Likelihood scores were graphed against generation time using Tracer v.l. 2.1. 
(Rambaut and Drummond 2005) to identify stationarity, and thus to determine how many generations must be 
discarded as bum-in, and whether or not more generations were required to be run. In order to insure that the 
analyses had found the optimal arrangements, they were performed twice for each data group and the station¬ 
arity levels were compared for convergence. When the different analyses reached stationarity and the topolo¬ 
gies were congruent, the resultant trees were combined using a majority-rule consensus tree in PAUP ver.4 
(Swofford 1998). Congruence for each branch indicated the posterior probability (PP). Using the criterion of 
a=5%, clades were considered to be significantly supported when PP =95% (Wilcox et al. 2002; Reeder 
2003). 

Choosing the outgroup 

In a preliminary analysis, the trees were rooted utilizing sequences of S. grammicus Wiegman and S. 
megalepidurus Smith, which are the first and second outgroups of the torquatus group (Wiens and Reeder 
1997), respectively. However, the furthermost external group contributed less in terms of character states and 
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rooting information, and introduced errors into the analysis (Lyons-Weiler, et al. 1998; Nylander 2001; Sand¬ 
erson and Shaffer 2002). Therefore, we added one more taxon to the first outgroup (5. heterolepis Boulenger), 
thus breaking the long branch leading to the external group and adding balance to the topology (Swofford and 
Olsen 1990; Smith 1994). 


Results 

Sequences 

Sequence data from 54 lizards belonging to 25 taxa were assembled. We could not amplify 12S from 5. 
bulleri 1 and 5. bulleri 2. We obtained 912 and 925 bp of the ribosomal genes 12S and 16S respectively, and 
709 bp of encoding ND4 plus 184 bp of the adjacent tRNAs. After alignment, a matrix of 2701 characters was 
obtained, of which 1789 were constant, 304 were variable but not phylogenetically informative, and 599 were 
potentially phylogenetically informative. 

Phylogenetic analyses 

Analyses of the separate genes typically resulted in congruent topologies. The relationships are similar for 
those obtained in the total evidence analyses, with the exception of two incongruent nodes, that were weakly 
supported (the trees are not shown, but they are available upon request). First, in the analyses of 16S and ND4 
(and also in the combined analysis), 5. insignis Webb was resolved as the sister taxon of a clade formed by 
(((5. spl + (S. t. melanogaster 1 + S. t. melanogaster 2)) + (S. t. torquatus Wiegmann + S. t. binocularis 
Dunn)) + ((5. bulleri 2 + 5. bulleri 3) + S. bulleri 1)), but with low BSPs (16s and ND4: BSP=57) and high 
and moderate PPs (16s: PP=100; ND4: PP=94). Alternatively, 12S recovered 5. insignis as the sister taxon of 
a clade formed by ((5. j. jarrovii 11a + 5. j. jarrovii lib) + ( S.jarrovii 10 + 5. lineolateralis Smith)) with a 
weak support (BSP<50, PP=58). Second, 5. ornatus caeruleus Smith was in a polytomy in the analyses of 16S 
and ND4, but 12S (like in the combined analysis) recovered it as the sister species of a clade formed by (((5. 
cyanogenys 1, 2, 96) + (5. plioporus Smith, 1+5. plioporus 2)) + (5. j. cyanostictus Axtell and Axtell, la + 5. 
j. cyanostictus lb)) with a weak and strong support (BSP<50, PP=98). 

The MP of the combined data found 30 most parsimonious trees (MPTs: length=2254, CTO.524, 
RI=0.733) and the strict consensus tree is shown in Figure 2. For the BI of the combined data, the first 2000 
generations were discarded as the burn-in. The strict consensus of the MPTs and the Bayesian majority-rule 
probability tree of 18001 trees were congruent in their relationships, although MP recovered 5. o. caeruleus, 
one population of 5. oberon Smith and Brown (5. oberon 21b) and one population of S. minor (S. minor 13a) 
as a polytomy. The better resolved BI tree is our preferred phylogenetic hypothesis and is presented in Figure 
3 along with PP and BSPs support above and below the branches, respectively. This hypothesis (Fig. 3) shows 
two strongly supported basal clades, A and B (BSP=100 and PP=100). Clade A includes the subspecies of 5. 
torquatus plus 5. bulleri, 5. insignis, 5. sp. 1, 5. linoelateralis and 5. jarrovii. whereas Clade B includes the 
remaining species. 

Clade A has two strongly supported subclades (BSP=100, PPM 00). In one subclade 5. lineolateralis was 
resolved within populations of S.jarrovii and this association was strongly supported (BSP=100, PPM 00). In 
the second subclade, 5. insignis was the sister taxon to all other species (BSP=52, PP=99). The next node of 
this subclade resolved the populations of 5. bulleri clade ((5. bulleri 1 + (5. bulleri 2 + 5. bulleri 3)) with 
strong support (BSPM00, PPM00). Sceloporus sp. 1 was resolved within the subspecies of 5. torquatus as 
((5. sp. 1,5. t. melanogaster ) + (5. t. torquatus, 5. t. binocularis )) with moderate and strong support (BSP=72, 
PPM00). Sceloporus torquatus subespecies were the sister group of 5. bulleri (BSP=95, PPM00). 

Within Clade B, Clade C was strongly supported (BSPM00, PP=95). Sceloporus mucronatus aureolus 
Smith was the sister of (5. m. omiltemanus Gunter + 5. macdougalli Smith and Bumzahem) and this clade 
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received moderate support (BSP=67, PP=95). Surprisingly, nominate S. m. mucronatus Cope was resolved 
with strong support (BSP=95, PP=100) in Clade G as the sister taxon of ( S . suguillatus Smith + S. poinsettii 
Baird and Girard). Therefore, the subspecies of S. mucronatus did not form a monophyletic group. 



S. cyanogenys 1 
S. cyanogenys 2 
S. cyanogenys 96 
S. s. plioporus 1 
S. s. plioporus 2 
S. cyanostictus 1 a 
S. cyanostictus 1 b 
S. oberon 24a 
S. oberon 28a 
S. oberon 23 
S. oberon 29a 
S. o. ornatus 
S. oberon 27b 
S. oberon 21b 
S. o. caeruleus 
S. minor 6a 
S. minor 4 b 
S. minor 5 
S. minor 3 
S. minor 8 
S. minor 1 
S. minor 1 5a 
S. minor 1 4a 
S. minor 17 
S. s. serrifer 1 
S. s. serrifer 2 
S. s. prezygus 1 
S. s. prezygus 2 
S. minor 13a 
S. poinsettii 
S. sugillatus 
S. m. mucronatus 
S. d. intermedius 1 
S. d. intermedius 2 
S. d. dugesii 
S. m. omiltemanus 
S. macgdougalli 
S. m. aureoius 
Sceloporus sp. (1) 
Sceloporus sp. (2) 
Sceloporus sp. (3) 

S. t. melanogaster 2 
S. t. melanogaster 1 
S. t. torquatus 
S. t. binocularis 
S. bulleri 3 
S. bulleri 2 
S. bulleri 1 
S. insignis 2 
S. insignis 1 
S. jarrovii 11 a 
S. jarrovii 11 b 
S. jarrovii 1 0 
S. lineolateralis 
S. grammicus 
S. heterolepis 


FIGURE 2. Srict consensus of 30 trees from the parsimony analysis based on 12S 16S and ND4 mtDNA sequences 
(length=2254, CI=0.524, RI=0.733). Bootstrap proportions > 50 % are indicated above the branches. 


SCELOPORUS TORQUATUS SPECIES-GROUP 


Zootaxa 1609 © 2007 Magnolia Press ■ 61 
































































Clade D contains Clade E as the sister group of Clade F. In Clade E, the monophyly of S. dugessi Bocourt 
is well supported (BSP=81, PP=100). This result differed from that of Wiens and Reeder (1997) who, with 
weak support, placed S. d. dugesii as a sister taxon of S. poinsettii. 



S. cyanogenys 1 
S. cyanogenys 2 
S. cyanogenys 96 
S. s. plioporus 1 
S. s. plioporus 2 
S. cyanostictus 1 a 
S. cyanostictus 1 b 
S. o. caeruleus 
S. oberon 24a 
S. oberon 28a 
S. oberon 23 
S. oberon 21b 
S. oberon 29a 
S. o. ornatus 
S. oberon 27b 
S. minor 6a 
S. minor 4b 
S. minor 5 
S. minor 3 
S. minor 8 
S. minor 1 
S. minor 15a 
S. minor 14a 
S. minor 17 
S. minor 13a 
S. s. serrifer 1 
S. s. serrifer 2 
S. s. prezygus 1 
S. s. prezygus 2 
S. poinsettii 
S. sugHiatus 
S. m. mucronatus 
S. d. intermedius 1 
S. d. intermedius 2 
S. d. dugesii 
S. m. omiltemanus 
S. macgdougalli 
S. m. aureolus 
Sceloporus sp. (1) 
Sceloporus sp. (2) 
Sceloporus sp. (3) 

S. f. melanogaster 2 
S. f. melanogaster 1 
S. t. torquatus 
S. t. binocularis 
S. bulleri 3 
S. bulleri 2 
S. bulleri 1 
S. insignis 2 
S. insignis 1 
S. jarrovii 11 a 
S. jarrovii 11 b 
S. jarrovii 1 0 
S. lineolateralis 
S. grammicus 
S. heterolepis 


FIGURE 3. Bayesian inference tree based on 12S 16S and ND4 mtDNA sequences. Posterior probabilities > 50% and 
boostrap proportions > 50 % (from the parsimony analysis) are indicated above and below the branches, respectively. 
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Clade F had two primary groups, clades G and H, and in turn, Clade FI contained clades I and J. Clade I 
(BSP=74, PP=88) consisted of two subspecies of S. serrifer Cope plus S. minor, S. s. serrifer was the sister to 
S. s. prezygus Smith (BSP=100, PP=100) and together they formed the sister group of S. minor (Fig. 3). How- 
ever, the subspecies of S. serrifer were not recovered as a monophyletic group because S. s. plioporus was 
resolved as the sister taxon of S. cyanogenys in Clade J (BSP<50, PP=100), an arrangement that agreed with 
the morphological analysis of Olson (1987). 

Incidences of non-monophyly occurred in Clade J. Sceloporus ornatus ornatus Baird, branched off from 
within S. oberon and S. o. caeruleus was the sister group of S. cyanostictus, S. s. plioporus, and S. cyanogenys. 
The phylogenetic relationships of S. cyanogenys, S. cyanostictus, S. oberon and S. minor are in discordance 
with the analysis of Wiens et al. (1999). We recovered moderate and strong support (BSP=71, PP=97) for the 
placement of S. oberon as sister taxon of S. cyanogenys, S. cyanostictus, and S. plioporus. In contrast, Wiens 
et al. (1999) reported a weakly supported subclade where S. minor and S. oberon were the sister group of S. 
cyanogenys and S. cyanostictus. 


TABLE 2. Data partitions, the best models of sequence evolution according to the BIC method and the number of char¬ 
acters of each partition used in the Bayesian inference analysis. 


Partition 

Model 

Number of characters in partition 

12S 

TrN+I+G 

912 

16S 

TrN+I+G 

925 

ND4 1 st codon 

HKY+G 

237 

ND4 2 nd codon 

TrN+I+G 

236 

ND4 3 rd codon 

TrN+G 

236 

tRNAs 

K81 uf+G 

184 


TABLE 3. Values of the parameters, estimated using the BIC method of Bayesian Inference for the different data groups. 



Substitution rates 





Ti/tv 

ratio 

Site rates 

Nucleotide frecuencies 


A<->C 

A<->G 

A<->T 

C<->G 

C<->T 

G<->T 


I 

A 

A 

C 

G 

T 

12S 

1.0000 

5.9702 

1.0000 

1.0000 

9.3163 

1.0000 

- 

0.3641 

0.4946 

0.3778 

0.2361 

0.1649 

0.2212 

16S 

1.0000 

2.7657 

1.0000 

1.0000 

6.1090 

1.0000 

- 

0.4981 

0.5931 

0.3735 

0.2414 

0.1658 

0.2194 

P'codon 

- 

- 

- 

- 

- 

- 

4.6456 

- 

0.3055 

0.3425 

0.2571 

0.1717 

0.2287 

2” d codon 

1.0000 

19.4007 

1.0000 

1.0000 

1.3595 

1.0000 

- 

0.5900 

0.3189 

0.1632 

0.2936 

0.1454 

0.3978 

3 rd codon 

1.0000 

32.0663 

1.0000 

1.0000 

10.6605 

1.0000 

- 

- 

2.5339 

0.4534 

0.3046 

0.0478 

0.1942 

tRNAs 

1.0000 

2.9227 

0.1462 

0.1462 

2.9227 

1.0000 

_ 

_ 

0.4739 

0.3709 

0.1958 

0.1575 

0.2758 


Discussion 

In this study, significantly and moderately supported relationships were obtained for all species and subspe¬ 
cies of the torquatus group. In general, the relationships agreed with those of Wiens et al. (1999), although 
some of the relationships recovered by Wiens and Reeder (1997) differed. Unlike Wiens and Reeder (1997), 
the monophyly of the subspecies of S. torquatus and their relationships with S. bulleri and S. insignis were 
strongly supported. S. sp. 1 from Jalisco was the sister taxon of S. t. melanogaster. Similarly, S. lineolateralis 
was resolved as the sister taxon of S.jarrovi, as suggested by Sites et al. (1992). 
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The non-monophyly of S. mucronatus was corroborated according to Wiens and Reeder (1997) and S. m. 
mucronatus was supported as being the sister taxon of ( S. sugillatus + S. poinsettii ) (BSP=86, PP=100). Sce- 
loporus m. aureolns was the sister taxon of (S. macdougalli + S. m. omiltemanus) (BSP=67, PP=95). In con¬ 
trast, Wiens and Reeder (1997) reported that S. macdougalli was the sister taxon of ( S. m. aureolus + S. m. 
omiltemanus), albeit with weak support. We found support for the monophyly of S. dugesii (BSP=81, 
PP=100), a relationship that also differs from the analysis of Wiens and Reeder (1997). They resolved S. d. 
dugesii as the sister taxon of S. poinsettii, and S. d. intermedins as the sister taxon to all other species of the 
torquatus group. 

In the study of Wiens and Reeder (1997), S. s. serrifer, S. s. prezygus, and S. cyanogenys, were not found 
to be sister taxa. However, Wiens and Reeder treated that result with caution, given the low branch support. 
Similarly, our hypothesis (Fig 3) did not resolve these taxa as being a monophyletic assemblage. This finding 
contrasts with the morphological evidence of Olson (1987), who proposed that S. cyanogenys was a subspe¬ 
cies of S. serrifer. 

Olson (1987) associated S. s. plioporus (not included by Wiens and Reeder 1997) with S. cyanogenys. In 
our study, S. serrifer plioporus was the sister taxon of S. cyanogenys and this association received strong sup¬ 
port (BSP=100, PP=100). This association also has geographical support. Whereas both S. s. serrifer and S. s. 
prezygus occur in southeastern Mexico, S. s. plioporus principally inhabits southern Tamaulipas and a small 
portion of northern Veracruz (Fig. 1). This is south of the distribution of S. cyanogenys. Olson’s results (1987) 
as well as ours show that S. s. plioporus forms the southern part of a morphological cline of S. cyanogenys, 
and should be considered as the same species. A single morphological characteristic typically differentiates S. 
s. plioporus from S. cyanogenys. In S. s. plioporus, the supraocular scales are complete and separated from the 
parietals by a row of intervening small scales. Alternatively, in S. cyanogenys, the supraocular scales are 
divided and in contact with the parietal scales. Within populations of S. s. plioporus in Tamaulipas, both mor¬ 
phological conditions exist. The percentage of individuals with divided supraocular scales increases north¬ 
wardly. Similarly, in some individuals, the supraoculars contact with the parietals, and in others they do not. 
The percentage of individuals that have supraoculars contacting the parietals diminishes northwardly. Unfor¬ 
tunately, we could not locate any population of S. serrifer from Veracruz (Smith 1939; Stuart 1970 and Olson, 
1987). A large percentage of Veracruz has suffered deforestation and been subjected to other types of ecologi¬ 
cal modification. For that reason, we cannot genetically determine whether these populations are more closely 
associated with S. cyanogenys or with S. serrifer of southeastern Mexico. 

Wiens and Reeder (1997) resolved the two subspecies of S. ornatus as sister taxa. However, no molecular 
data were available for S. o. ornatus and the association was weakly supported. In contrast, the two subspecies 
were not recovered as sister taxa in our study. Sceloporus o. ornatus occurs in Coahuila, and is geographically 
close to populations of S. oberon (Fig. 1). Although possible, we do not believe that our results are the conse¬ 
quence of a recent invasion or introgression of the maternal genotype of S. oberon into S. o. ornatus. If migra¬ 
tion was involved, then we would expect S. o. ornatus to be more closely related to the geographically closest 
population, that of S. oberon 27 from Coahuila (see Figure 1 and Table 1). However, S. o. ornatus appeared as 
the sister group of the geographically furthermost population from Nuevo Leon ( S. oberon 29). Moreover, 
these taxa occur in very different environments. Whereas S. oberon occurs in oak woodlands, S. o. ornatus 
lives at lower altitude in desert regions. While S. oberon exhibits dark colors on its back, S. o. ornatus is yel¬ 
low and light blue. With respect to the dorsal scales, S. oberon has relatively large scales, averaging 37.5 
around the body, but S. o. ornatus averages 55, smaller scales. Our tree leaves three possible options to con¬ 
sider: 1) S. oberon and S. o. ornatus form a single species; 2) S. oberon contains at least three cryptic species; 
or 3) the non-monophyly owes to incomplete lineage sorting. An evaluation of highly variable nuclear genes 
could differentiate between these possibilities. However, for the time being, we prefer the first option and con¬ 
sider the taxa to be conspecific. The phylogenetic relationships of S. o. caeruleus are still not satisfactorily 
resolved, and we believe that more detailed studies are necessary. 
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The extensive variation in coloration between individuals in the torquatus group may reflect sexual selec¬ 
tion (Wiens et al. 1999). Regardless, environmental characteristics might also play a very important role, par¬ 
ticularly in the number and size of the dorsal scales, given that scales are involved in thermoregulation and 
humidity exchange (Soule and Kerfoot 1972; Fox 1975). 

Taxonomy of the torquatus group 

In order to obtain a taxonomy that reflects phylogenetic history, a number of taxonomic changes are nec¬ 
essary. The following modifications are proposed: 

1) Sceloporus mucronatus should be treated as a monotypic species. The subspecies S. mucronatus mucr- 
onatus should not be recognized. 

2) The subspecies Sceloporus mucronatus aureolus should be elevated to full species status as Sceloporus 
aureolus [new combination]. 

3) The subspecies Sceloporus m. omiltemanus should be elevated to full species status as S. omiltemanus 
[new combination]. 

In this study, we showed molecular evidence for the non-monophyly of S. mucronatus subspecies, which 
indicates a discordance between morphological and mtDNA species limits. The main differences between S. 
mucronatus subspecies have traditionally been identified as some patterns on the coloration, the number of 
dorsal scales and femoral pores (Smith 1939). In S. m. mucronatus dorsal scales are 27 to 30 with 11 to 17 
femoral pores on each side; in S. m. omiltemanus dorsal scales are 30 to 38 with 12 to 16 femoral pores, and in 
S. m. aureolus dorsal scales are 32 to 36 with 12 to 16 femoral pores. Nevertheless, due to wide moiphological 
overlapping between species, no consistent diagnostic characters have been observed. 

4) Sceloporus oberon should be synonymized into Sceloporus ornatus. Sceloporus ornatus Baird, 1859 
has priority over S. oberon (S . jarrovii oberon Smith and Brown, 1941). Although recognition of subspecies 
has become controversial, S. ornatus ornatus could continue to be recognized. If so, then populations pres¬ 
ently known as S. oberon should be referred to as S. ornatus oberon [new combination]. We recognize that this 
arrangement results in a paraphyletic taxonomy for the subspecies. 

5) Sceloporus ornatus caeruleus should be elevated to full species status as S. caeruleus [new combina¬ 
tion]. 

As in S. mucronatus, we observed discordances between morphology and mtDNA data. According to the 
molecular phylogeny of Wiens et al. (1999) S. jarrovii oberon and some northern populations of S. j. minor 
are synonymized in S. oberon, despite the differences in coloration of these two taxa. Wiens et al. (1999) sug¬ 
gested that the differences in dorsal coloration in the populations of S. oberon may reflect sexual selection. 
Furthermore, in our study we also found that S. o. ornatus and S. oberon conforms an evolutionary species, 
despite the differences in coloration and scutelation. The populations of S. oberon have between 34 to 46 dor¬ 
sal scales, whereas S. o. ornatus have between 55 to 63 dorsal scales with a complex coloration pattern (Smith 
1939). The differences in the number of dorsal scales may be due to habitat, as was pointed out in a previous 
paragraph. Habitat influence may explain the morphological similarities between S. o. ornatus and S. o. caer¬ 
uleus which has a high number of dorsal scales (47 to 53) and also occurs in semi-desert habitats, but without 
a close phylogenetic relationship. 

6) Sceloporus serrifer plioporus Smith, 1939 from southern Tamaulipas, should be synonymized into S. 
cyanogenys Cope, 1885. The taxonomic status of populations in Veracruz remains uncertain. The original 
morphological difference between putative populations of S. s. plioporus and S. cyanogenys, was the divided 
supraoculars scales in the latter (Smith 1939). However, on closer inspection, these differences are not sup¬ 
ported (Olson, 1987) because the percentage of individuals with divided supraoculars scales increases north¬ 
wardly. 

7) Sceloporus dugesii should be recognized as being monotypic, instead of having two subspecies S. d. 
dugesii and S. d. intermedius. Despite Wiens and Reeder (1997) found some weakly supported morphological 
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differences between S. d. dugesii and S. d. intermedins. The main diagnostic character between these two taxa, 
the presence of head scales microscopically rugose in S. d. dugesii (Smith, 1939), is not a fixed character (see 
morphological matrix in the study of Wienes and Reeder 1997), and it could be chosen on a small sample size 
basis. 

8) Sceloporus lineolateralis Smith, 1936 should be synonymized into Sceloporus jarrovii Cope, in Yar¬ 
row, 1875, but potentially recognized as the subspecies S. jarrovii lineolateralis [new combination]. Unfortu¬ 
nately S. j. jarrovii lacks fixed diagnostic morphological characters (Wiens and Penkrot 2002). The characters 
early identified by Smith (1939) like diagnostic of S. j. jarrovii (e. g., the first canthal seldom forced above 
canthal ridge by contact of second canthal and subnasal, prefrontals in contact, color pattern etc.), exhibit 
some intraespecific variation even in other populations. Some authors have similarly reported that S. lineolat¬ 
eralis and S. j. jarrovii intergrade with each other based on morphological characters (Webb and Hensley 
1959; Chrapliwy 1964; Wiens et al. 1999). The previous studies along with our molecular results indicate the 
conspecificity between S. lineolateralis and S. j. jarrovii. 
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